Stochastic Weighted Fractal Networks 



O 

(N 

PL, 



Timoteo Carletti 
Departement de Mathematique, 
Facultes Universitaires Notre Dame de la Paix 
8 rempart de la vierge B5000 Namur, Belgium 
timoteo. carletti ©fundp. ac. be 
(Dated: February 1, 2010) 

In this paper we introduce new models ol complex weighted networks sharing several properties 
with fractal sets: the deterministic non-homogeneous weighted fractal networks and the stochastic 
weighted fractal networks. Networks of both classes can be completely analytically characterized in 
terms of the involved parameters. The proposed algorithms improve and extend the framework of 
weighted fractal networks recently proposed in [3!|. 

PACS numbers: 89. 75. He Complex networks, 05.45.Df Fractals. 05. 40. -a Stochastic processes 
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I. INTRODUCTION 

Fractal structures are ubiquitous in nature, coast- 
lines river networks H y], snowflakes Q, grow- 
ing colonies of bacteria Qa TTL m ammalian lungs |cj-ll2|. 
mammalian bloody vessels |12| . just to mention few of 
them • But also mankind artifacts can exhibit fractal 
features, for instance fractal antenna fl5j or fluctuations 
in markets prices [lt| . 

A distinction can be made between mathematical or de- 
terministic fractals [TtJ for which a complete geometric 
description can be provided using simple tools such as ho- 
motheties, rotations and copying, and random or pseudo 
fractals [13| found in nature, being the latter character- 
ized by exhibiting fractal properties, for instance self- 
similarity, only when statistical averages are computed, 
because unavoidable fluctuations and errors can alter the 
regular-geometric patterns. Moreover such scale invari- 
ance should be limited to a finite range of scale lengths 
because of physical constraints. 

It is worth remarking that some of these physical frac- 
tals have functionalities, e.g. transportation of gases 
in mammalian lungs, or charges in fractal antenna, one 
can thus improve the geometrical description by includ- 
ing flows and growths constraints. Networks are there- 
fore the most natural and useful tool to describe such 
growing complex structures with flows constraints. We 
thus hereby propose the Stochastic Weighted Fractal Net- 
works, SWFN for short, a new class of complex networks 
whose construction is directly inspired by such physical 
fractal structures. 

Star ting from the pioneering works of Erdos and 
Renyi [lj|, network theory is nowadays a research field 
in its own (l9l . [20| and the scientific activity is mainly 
devoted to construct and characterize complex networks 
exhibiting some of the remarkable properties of real net- 
works, scale-fr ee l2ll. small-world j22j . communities [23j . 
weighted links [24T - l29j . just to mention few of them. 

In recent years we observed an increasing number of 
papers were authors proposed models of deterministic 
(pseudo) fractal networks (30rl39j exhibiting scale-free 
and hierarchical structures. In a limited number of cases, 



models presented also a stochastic component |4CH43| . 

The aim of the SWFN hereby introduced, is to provide 
a framework that could be used to (re) analyze flows on 
natural fractal structures using standard tools of trans- 
port theory on networks. Moreover SWFN share with 
physical fractals several interesting properties, for in- 
stance the self-similarity or the self- affinity, the presence 
of hierarchical structures and a stochastic growth pro- 
cess. Actually this allows us to generalize in a unifying 
scheme some of the above mentioned models existing in 
the literature. 

The SWFN are constructed via a stochastic process 
and we are thus able to analytically characterize their 
topology as a function of the parameters involved in 
the construction, using expectations obtained construct- 
ing several replicas. 

Let us conclude this introduction with two remarks. 
First of all we named our models "fractal" networks in- 
stead of "pseudo fractal" , because some of the topologi- 
cal properties of SWFN depend on the fractal dimension 
of some underlying fractal set, whose value ranges all the 
positive real numbers, without any limitation. Second we 
rather prefer to talk about "stochastic" networks to em- 
phasize the stochastic growth process instead of the ran- 
domness of some topological quantities; let us also stress 
that in the network theory " randomness" has a precise 
meaning that we cannot directly apply to this case. 

The paper is organized as follows. In the next section 
we will introduce and study a deterministic model, that 
generalize the one proposed in [39j . and that will serve 
as the basic building block to construct the SWFN in 
Section [TTT1 Then we conclude with some possible appli- 
cations we sum up and draw our conclusions 



II. DETERMINISTIC WEIGHTED FRACTAL 
NETWORK 



According to Mandelbrot [l7j "a fractal is by definition 
a set for which the Hausdorff dimension strictly exceeds 
the topological dimension". One of the most amazing 
and interesting feature of fractals is their self- similarity 
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FIG. 1: The map T s ,t, a - On the left a generic initial graph G 
with its attaching node a (red on-line) and a generic weighted 
edge w £ G (blue on-line). On the right the new graph G' 
obtained as follows: Let G' 1 ' , . . . , G' s ' be s copies of G, whose 
weighted edges (blue on-line) have been scaled respectively by 
a factor /i, . . . ,/ s , and let us denote by for i = 1, . . . , s, 
the node in G' % ' image of the labeled node a £ G, then link 
all those labeled nodes to a £ G through edges of unitary 
weight. The connected network obtained in this way will be 
by definition the image of G through the map: G 1 = T s ,f, a {G). 



or self -affinity [44|,|45|], namely looking at all scales we can 
find conformal or stretched copies of the whole set; this is 
actually the idea used to build up fractals as fixed point 
of Iterated Function Systems [HHIt}, IFS for short. Such 
fractals have a Hausdorff dimension completely charac- 
terized by the number of copies and the scaling factors of 
the IFS. Let us observe that in this case this dimension 
coincides with the so called similarity dimension [47j | . 

Recently, author proposed [39j a new general frame- 
work aiming to construct weighted networks with some 
a priori prescribed topology depending on the two main 
parameters: the number of copies and the scaling fac- 
tors, hence on the fractal dimension of the "underly- 
ing" IFS fractal. The aim of this Section is to generalize 
such construction to obtain a larger class of networks; 
moreover exploiting the iterative construction we will be 
able to completely and analytically describe the network 
topology in terms of node strength distribution, average 
(weighted) shortest path and (weighted) clustering coef- 
ficient. 

Let us fix a positive integer s > 1 and s real num- 
bers /i,...,/ s E (0,1) and let us consider a (possibly) 
weighted network G composed by TV nodes, one of which 
has been labeled attaching node and denoted by a. We 
then introduce a map, Ts,f, a , depending on the parame- 
ters s, f = (/i, . . . , f s ) and on the labeled node a, whose 
action on networks is described in Fig. [TJ 

So starting with a given initial network Go we can con- 
struct a family of weighted networks (Gk)k>o iteratively 
applying the map 7^,f, : G k = T s ,t, a {G k -i). 

This construction improves the one recently pro- 
posed [3!| by avoiding the introduction of an extra node, 
moreover it offers a unifying framework where several 
constructions presented in literature can be included 
and generalized, e.g. the model presented in [32[ with 
m = 3 can be mapped into to the WFN with s = 3, 



f = (1,1,1), i.e. no weights, and Go = Finally this 
deterministic construction will be the basic brick to de- 
velop the stochastic network introduced in the following 
Section Hill 

Given Go and the map 7~ s .f,a we are able to completely 
characterize the topology of each Gk for k > 1 and also 
of the limit network Goo , defined as the fixed point of the 
map, T s ,t,a{G oo) = Goo. 



A. Results 

The aim of this section is to describe the topology of 
the graphs Gk for all k > 1 and Goo, by analytically 
studying their properties such as the average degree, the 
node strength distribution, the average (weighted) short- 
est path and the (weighted) clustering coefficient. 

At each iteration step the graph Gk grows as the num- 
ber of its nodes increases according to 



N k = (s + l) k N , 



(1) 



being TVo the number of nodes in the initial graph, while 
the number of edges satisfies 



Ek = (s + l) k (E + l)-l 



(2) 



being Eq the number of edges in Go- Hence in the limit 
of large k the average degree is asymptotically given by 



Ek _ 
N k k- 
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(3) 



Let us denote the weighted degree of node i £ Gk, also 
called node strength pTjj . by = ^2j w ij ■> being 
the weight of the edge (ij) £ Gk', then using the recursive 
construction, we can explicitly compute the total node 



strength, W k = J2i w 
"2s 



(k) 



F 



and easily show that 
((F + l) k - 1) + (F + l) k W 



(4) 



being F = X^=i fj- Let us observe that using the hy- 
pothesis fj < 1, it trivially follows that F < s, hence we 
can conclude that the average node strength goes to zero 
as k increases: Wk/N^ — > 0. 



B. Node strength distribution. 

Let gk (x) denote the number of nodes in Gk that have 
strength cj| = x and let us assume go to have values in 
some finite discrete subset of the positive reals, namely: 



go(x) > if and only if x 6 {x\, . 



, Xj 



otherwise go(x) = 0. Using the property of the map 
7s, f, a we get that after k steps of the construction the 
nodes strengths have been rescaled by a factor f kl . . . f ka , 
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FIG. 2: Node Strengths Distribution. Plot of the renor- 
malized node strengths distribution D~ 1 log 10 gk{x), where 
D — — s log s/ log(/i ...f a ). Symbols refer to : □ the finite 
approximation Gn with 3145728 nodes of the WFN with 

s = 3, f = (1/^2,1/^3,1/^5) and G =<J; O the finite 
approximation Gg with 3359232 nodes of the WFN s = 5, 
f = (1/VE, 1/^,1/^3,1/^7,1/713) and Go =•-•. The 
reference line has slope —1, linear best fits (data not shown) 
provides a slope -1.037 ± 0.04 and R 2 = 0.798 for □ and 
-1.00 ± 0.03 and R 2 = 0.8382 for Q. 



where the non-negative integers ki do satisfy k\ + ■ • ■ + 
k s < k. Because this can be done in k\/(ki \ . . . k s \) pos- 
sible different ways, wc get the following relation for the 
node strength distribution for the network Gk- 



</*•(/• ' •••./: •<•; 



fc! 



fei!...fc s ! 



go(x) with ki + • • • + k s < k . 

(5) 

After sufficiently many steps and assuming that the main 
contribution arises from the choice k\ ~ • • • ~ k s ~ k/s, 
we can use Stirling formula to get the approximate dis- 
tribution (see Fig. [2]) 



logg k (x) 
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log(/i.../ s ) 
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C. Average weighted shortest path. 

By definition the average weighted shortest path [2(| of 
the graph Gk is 



where 



A fc = 



A, 



N k (N k - 1) 



(fc) 



(7) 



(8) 



being the weighted shortest path linking nodes i and 
j in Gk- Taking advantage of the recursive construction 
and adapting the ideas used in [Hj], we get the following 
recursive relation for A^ 

A fe = (F+ l)A k -i + 2s(F + l^fc-iA^T + 2s 2 AT fe 2 _ 1 , 

(9) 

i.e. the sum of 



where wc introduced A 



(a fc ) 



'-fc SieGfc Pia k i 

all weighted shortest paths ending at the attaching node, 
afc G Gk- We can prove that for large k the asymptotic 
behavior of A[, is given by 



A 



(afc) 



fc— >oo S 



(10) 



and thus the recursive relation ^ can be explicitely 
solved to provide the following asymptotic behavior in 
the limit of large k (see Fig. [3]) 



A* 



2s 2 (s + l) 



(11) 



fc^oo (s-F)[(l + s) 2 - (1 + F)] ■ 

One can explicitly compute the average shortest path, 
Ik-, formally obtained by setting fi = ■ ■ ■ = f s = 1 in the 
previous formulas ([7]) and ([S]). Hence slightly modifying 
the results presented above we can prove that asymptot- 
ically we have (sec Fig. 0]) 



2.s 



fc^oo (1 + s) log(s + 1) 



los 



Nk 
N 



(12) 



where growth law of Nk given by ([1} has been used. Thus 
the network grows unbounded with the logarithm of the 
network size, while the weighted shortest distances stay 
bounded. 



D. Clustering coefficient. 



so the nodes strength distribution follows a power law. 
Let us observe that in the case of homogeneous scaling 
i.e. all fj equal to some / G (0, 1), one can prove [39| 
that Eq. © reduces to log<7fc(a;) ~ — df ract log a; where 
dfract = - logs/ log / is the fractal dimension of the un- 
derlying IFS fractal. 



The clustering coefficient [13, [22j of the graph Gk is 

defined as the average over the whole set of nodes of the 

(fc) 

local clustering coefficient c- , namely < Ck >= Ck/Nk, 

where Ck=T,ieG k c i k) ■ 

Because of the construction algorithm new triangles 

are created in the network "boundary" while their number 



4 



o.s- 



■0.6- 



0,4- 



0.2 - 



9- 
8- 
7- 
6- 
5- 
4- 
3- 
2- 
1'A 




10 



10 



10' 
N k /N 



10 



10 



FIG. 3: The average weighted shortest path. Plot of the 
renormalized average weighted shortest path Afc versus the 



iteration number k, where Afc = Aa 
F — fi + ---+f s . Symbols refer to 



(s-F)[(l+s)Ml+F)] 



and 



2s 7 (s+l) 

□ the WFN s = 3, 



f = (1/V2, 1/V3, 1/V5) and Go =<J; O the WFN s =5, 
f = (l/\/5, 1/VTT, 1/V3, 1/V7, l/\/l3) and G = the 
WFN s = 2, f = (1^/3,1/^/5) and G =<*; * the WFN 

s = 2, f = (1/V3, l/\/5) and Go 



doesn't change in the inner core, hence the local cluster- 
ing coefficient, at each step increases just by a factor s+1; 
thus after fc-interactions we will have Ck = (1 + s) fc Co, 

being Co = X^eGo c i°^ ^ ne sum °^ l° ca l clustering coef- 
ficients in the initial graph. We can thus conclude that 
the clustering coefficient of the graph is asymptotically 
given by: 



< c k > — > — . 



^1 



(13) 



On the other hand, one can introduce the links val- 
ues to weigh the clustering coefficient [HI, generalizing 
the previous relation, we can easily prove that weighted 
clustering coefficient of the graph is asymptotically given 
by: 



< 7fe > = 



I + F\ K C 



1 + s ) N fc->oo N, 



l-d 



(14) 



where d = l ? s ^}, + . F ) , that results smaller than one because 

log(l+s) ' 

of the assumption fj < 1. 



FIG. 4: The average shortest path as a function of the 
network size (semilog plot). Semilog plot of the renormalized 
average shortest path Ik versus the network size N^, where 
4 = 4 (3+1) ' 2 ° S(3+1) - Symbols are the same of Fig. H The 
reference line has slope 1. Linear best fits (data not shown) 
provides a slope 0.970±0.017 and R 2 = 0.9999 for □, 0.9654± 
0.05 and R 2 = 0.9997 for O, 0.97 ±0.02 and R 2 = 0.9998 for 
and 1.01 ± 0.03 and R 2 = 0.9993 for *. 



III. STOCHASTIC WEIGHTED FRACTAL 
NETWORKS 

The aim of this section is to present a class of complex 
weighted networks that grow according to a stochastic 
process and exhibit self-similar or self-affine structures, 
hereby named Stochastic Weighted Fractal Networks, for 
short SWFN, whose construction is directly inspired by 
the stochastic growth phenomena present in nature. The 
idea is thus to mimic the growth of fractal structures in 
nature where "possible errors" could modify regular pat- 
terns. 

So let us hypothesize that the growth process is the 
result of a stochastic process that selects the actual real- 
ization, i.e. the number of copies, between a number of 
different possibilities. Thus at each iteration the number 
of copies, s, is a stochastic variable distributed accord- 
ing to some probability distribution function p(s). Once 
the numerical value for s has been set, s real numbers 
fi , . . . , f s are drawn according to some probability dis- 
tribution function q(f) with values in (0,1). Finally a 



new network is constructed by applying 7^ (/, 
the actual network: 



,f.),a 



to 



G 



p{s) 



G (s) = T s 



ih,-,h): 



rG. 



(15) 



Remark. In the following we will assume the simpli- 
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fying working hypothesis that f\ = ■ ■ ■ = f s = tt/s, i.e. 
q(f) = S(f — a/s), for some given and fixed a £ (0, 1), 
but of course the model applies to more general cases. 

One can repeat the construction k times and thus ob- 
tain with probability p(sk) ■ ■ -p(si), starting from a net- 
work Go, a new network, denoted by G^ Sk ''"' S1 ': 

G(Sfc ""' Sl) = r Sfc ,(/ 1 (fc, ,..,/i: , ),a ° • ' " ° T su(f^,-Aa G ° ■ 

(16) 

The network growth results thus a stochastic process, 
hence we will describe the main topological network mea- 
sures in terms of expectations obtained repeating several 
times the construction. Of course we could also consider 
and compute higher order momenta, but the computa- 
tions become rapidly cumbersome, and thus we will non 
present these results except for some simple cases, such 
as the number of nodes. 



A. Results: SWFN 

At each step the number of nodes increases with re- 
spect the present ones, and the exact amount depends 
on the number of branches drawn. Starting from a net- 
work containing Nq nodes we get a new network with 
]\[( s i) = (l + si)Nq nodes with probability p(si). Iter- 
ating the construction, after k steps we can obtain with 
probability p(s k ) . . . p(si) a network with i\K Sfc '"" Sl ) = 
(1 + Sfc) . . . (1 + si)Nq nodes. Hence the expected value 
for the number of nodes in a network build after k itera- 
tions, is given by: 

<N k >= P(s k )--.p(si)N^'-^ (17) 

Sfc,...,si 

= 5>( Sfe )(l + Sfc ) J2 P{s k -i)...p{ Sl )N^-^-<^ 

Sfc SJt_l,...,Sl 

= (1+ < s >) < JV fc _i >, 

where we denoted by < s >= J2s k p( s k)s k the aver- 
age number of branches. We can thus conclude that the 
expected number of nodes increases exponentially 

< N k >= (1+ < s >) k N . (18) 

Using similar ideas one can prove that the variance of the 



number of nodes increases according to: 

a 2 Nk = N 2 [(l+<s>) 2 + a 2 ] k -(l+< S >) 2k N , (19) 

where a 2 is the variance of the distribution of number of 
branches. 

On the other hand the number of edges can increase, 

with probability p(s k ), in one iteration by E^ Sk S1 > = 

(1 + s k )E( Sk - 1 '---' Sl > + s k and thus the expected number 
of edges do satisfy 

<E k >=(l+<s>) k (E + l)-l. (20) 
These findings are exact in the case of infinitely many 
replicas, nevertheless numerical simulations presented in 
Fig. [5] and in Fig. [()] show the good agreement also for 
finitely many repetitions. 

Remark. The numerical simulations presented in the 
following will be obtained assuming for the branch number 
a Poisson distribution translated by one, more precisely 
to avoid a non zero probability of drawing zero branches, 
we drawn with probability p(k) = A k e~ x jk\ a non nega- 
tive integer k, and then we set the number of branches to 
s = k + 1, in this way we will get < s >= A + 1. a 2 = A 
and s > 1. 

Of course our findings are more general and do not rely 
on the particular choice for p(s). 

In a similar way we can compute the expected average 
degree after k steps, < (E/N) k >, and the expected av- 
erage node strength after k steps, < (W/N) k >, where 
W is the total node strength for the given network real- 
ization, to get (see Fig. [5]): 

(21) 

As we did in the previous section, we are able to analyt- 
ically study other relevant quantities such as the expected 
value for the weighted shortest path < X k >, defined for 
each network realization by ([?])• More precisely, starting 
from a network Go and applying iteratively the above 
construction we end up after k iterations with probabil- 
ity p(s k ) . . .p(si) to a network G < - Sfc, " ,Sl \ we can thus 
define the weighted shortest path for the given network 
realization by A < - Sfc ' ' :Sl - 1 = , A ' fc ' — \ > . ^ . Then using the 

recursive construction we get: 



A(sfe ,... >si) = (Fk + l)A B *-i Sl + 2s 2 (jy('*-i.--«i)) | 2sk {F k + l)A^- 1, "" ,8l ^ (a '- 1 -- ai) 

(1 + s fc ) 2 (A^-i<-< s i)) 2 (1 + s fe ) 2 (iVK-i s i)) 2 



F k + 1 

(i + s k y 



\(,s k -i,-,si) + 2 



Sk 



s k + 1 



2s fc (F fc + l) Vsfc si) 

(1 + Sk) 2 



(22) 



r 



where we used the growth rate of the num- bcr of nodes and ^ and we introduced 
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FIG. 5: Expected values for number of Nodes and number of Edges. Renormalized quantities : A 1 log(< Nk > /No) and 
A -1 log < E k > where A = log(l+ < s >). Symbols refer to : O the SWFN with parameters A = 4, a — 0.5 and Go = 

□ the SWFN with parameters A = 2, a = 0.5 and Go =<J; the SWFN with parameters A = 3, a = 0.8 and Go =<J. 

Expectations are obtained over 100 replicas. Left panel, the reference line has slope 1, linear best fits (data not shown) give 
0.9998 ± 0.03 R 2 = 0.9991 for O and 1.017 ± 0.008 R 2 = 0.9999 for □, 1.008 ± 0.005 R 2 = 1.000 for <>. Right panel, the 
reference line has slope 1, linear best fits (data not shown) give 0.9569 ± 0.06 R 2 = 0.9988 for O, 1-019 ± 0.03 R 2 = 0.9997 for 

□ and 1.06 ± 0.07 R 2 = 0.9955 for 0. 
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FIG. 6: Expected values for the average degree and the average node strength. Renormalized quantities : < (E/N)k > /D where 
D = (Eo + 1)/Nq. Symbols are the same of Fig. [5] Expectations made over 100 replicas. 

^(sfc,...,si) _ ^(sfc.---»si) /^Y(s fe ,...,si)_ average weighted shortest path satisfies the recurrence 

One can finally prove that the expected value for the equation: 



(A fe ) = (A fc _!) 




(23) 
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where we defined 
1 + F 




? p(fc) (TT^ 



Sk 



s k + 1 



and 



t-r 1 (l + s k y 



Under the simplifying assumption fi = ■ ■ ■ = f Sk = 
a/ Sk we get = a and thus we can simplify the previous 
equations and obtain (see Fig. [7]): 



< A fe > -> 

k— yoo 



s+lj I 1 - (a + f ) (1/(1 + s) 2 } \ (1 + s) 2 



2(a + l) 



1 



/ 1 - (a + 1) (1/(1 + s) 2 } 1 - (a + 1) (1/(1 + a)) 
(25} 
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FIG. 7: Expected values for the average weighted shortest path. 
Renormalized quantities: < >= L^ 1 < >, where L is 
the right hand side of Eq. (|25[) . Symbols are the same of 
Fig. [5] Expectations made over 20 replicas. 



One can consider the expected shortest path by formally 
set all the scaling factors equal to 1 and similar technics 
allow to conclude that (see Fig. [8]) 



FIG. 8: Expected values for the average shortest path as a 
function of the network size (semilog plot). Semilog plot of the 
renormalized expected average shortest path < £k > — < £k > 
jM versus the network size iVfc, where M is the right hand 
side of Eq. (|26p . Symbols are the same of Fig. [5] Expectations 
made over 20 replicas. Fhe reference line has slope 1, linear 
best fits (data not shown) provides 0.95 ± 0.01 with R 2 = 
0.9999 for □, 0.95 ± 0.07 R 2 = 0.9968 for O and 0.98 ± 0.01 
R 2 = 0.9999 for 0. 



<4 > 




1 



1 - <l/(s + 1)) log(l+ < a >) 

(26) 

Remark. Let us observe that in the case where only 
one value of s is possible, i.e. the probability distribution 
of the number of branches reduces to a 8 -distribution, 
p(s) = 8 St3 i, then the above result coincide with the ones 
presented for the WFN in Section [771 



log- 



IV. CONCLUSIONS 



No 



In this paper we proposed a unifying general frame- 
work for complex weighted networks sharing several 
properties with fractal sets, the Stochastic Weighted 
Fractal Networks. This theory, that generalizes to net- 
works the construction of physical fractals, allows us 
to build complex networks with a prescribed topology, 
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whose main quantities can be analytically predicted in 
terms of expectations and have been shown to depend 
on the fractal dimension of some underlying fractal; for 
instance the networks are scale-free, the exponent being 
the related to the fractal dimension of the underlying IFS. 
Moreover the SWFN share with fractals, the self-similar 
or self-affine structure. 

These networks exhibit the small-world property. In 
fact the average shortest path increases logarithmically 
with the system size; hence it is as small as the average 
shortest path of a random network with the same number 



of nodes and same average degree. On the other hand 
the clustering coefficient is asymptotically constant, thus 
larger than the clustering coefficient of a random network 
that shrinks to zero as the system size increases. 

As already observed [3!| the self-similarity property of 
the SWFN make them suitable to model real problems 
involving some kind of diffusion over the network coupled 
with local looses of flow, here modeled via the parameters 
/ < 1. Moreover the stochastic growth process allows us 
to introduce more realism in the construction. 
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